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Abstract 

Using a many-body theory, we calculate the change in differential conductance AdI /dV 
after adding an impurity either on a clean surface or inside an elliptic quantum corral. 
Using the same set of parameteres for both cases, the qualitative features of the voltage 
dependence of AdI/dV observed in recent experiments, are reproduced. 

Introduction In recent experiments a Co atom (acting as a magnetic impurity) was placed 
at one focus of an elliptic corral, and a depression in the differential conductance dl /dV as a 
function of voltage (a signature of the Kondo effect) was observed not only at this focus, but 
also at the other one jl|. The space dependence of AdI /dV ( dl/dV after substracting the 
corresponding result for the empty corral) reflects mainly the density of the eigenstate of two- 
dimensional free electrons confined in the corral which lies at the Fermi energy. This space 
dependence has been reproduced by several theories. The voltage dependence is more subtle 
and requires in principle a many-body calculation of the Kondo resonance. However, most 
theories treat the interactions in a phenomenological way. Many-body effects were included 
either by perturbation theory M or by numerical diagonalization in a restricted subspace 
H . Unfortunately, since the separation of the relevant eigenstates of the corral (those with a 
sizeable hybridization with the impurity) is larger than the Kondo temperature these 
numerical results cannot reproduce the observed line shape. 

Here we extend our previous approach to the clean surface and take into account the 
direct tunneling between the tip and the impurity in order to compare the line shape of 
AdI /dV for an impurity inside the corral or on a clean surface. 

The model The simplest description of the problem includes a non-degenerate highly cor- 
related impurity state and the surface states near the Fermi energy ef- For the clean Cu(lll) 
surface, the latter are uncoupled from bulk states for wave vectors inside the neck of the Fermi 
surface 0. Band structure calculations indicate that these states have a weight larger than 
60% in the first layer. The remaining weight is mainly due to other surface states at differ- 
ent energies Q. Taking a standard r~ 7 ^ 2 distance dependence between sp and d electrons, 
and considering that the impurity atom enters in the three-fold coordinated position at the 
surface, one obtains that the hybridization of the impurity with the surface layer is nearly 
17 times larger than that with the second layer. Since the density of sp states at the first 
layer is also larger Q , this information suggests that the direct hybridization of the impurity 
with bulk states can be neglected. Taking into account that the distance between Cu atoms 
is smaller than the Fermi wave length 2tv/Icf ~ 30A, the hybridization of the impurity with 
a surface state j can be approximated as VXtpj (Ri), where Ri is a two-dimensional vector 
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denoting the position of the impurity, <pj(r) is the wave function of state j normalized as 
J \ipj(r)\ 2 d 2 r = 1, A is a length scale which we take as the square root of the surface per Cu 
atom (A = 2.38 A), and V is an energy which would represent the hopping between one atom 
of the surface and the impurity in a tight-binding description. For the clean surface, j labels 
the allowed wave vectors of extended Bloch states in a large area with periodic boundary 
conditions, while for the mirage experiment, the <Pj{r) are localized inside the corral. 
The above considerations lead to the following Anderson model: 

H = + E d d\d a + Ud\d^d\di 

ja a 

+v\^y P] {R i )dU ] . (i) 

Here cj^ and dj. create an electron on the j th conduction eigenstate and the impurity respec- 
tively. 



Approximations and relevant equations At sufficiently low temperatures, the differen- 
tial conductance dl jdV when the probing tip is at position r, is proportional to the density 
of the mixed state f a (r) = \ Y] . ipj(r)cj a + qd a . 

^jp- ~ Pf (e F + eV) = - -ImG f (e F + eV) , (2) 

where G/(u) = ({fa-', fl))^ is the Green function of fa{r) and q is the ratio between tun- 
neling matrix elements between impurity and tip and between impurity and surface. It is 
different from zero only for r very near Ri. The change in dl /dV after adding the impurity 
is proportional to the corresponding change in pf. From the equations of motion: 

AG/(r,u) = (VG° c (r,R i ,w)+q)(VG°(Ri,r,Lj) + q)G d (w), (3) 

where G<j(w) = {(da-', dl))u is the impurity Green function and G° c is the conduction electron 
Green function for V — 0: 

G c (r,r ,uj )=> j— , 4 

^ — ' LJ + 10 — Ej 

j 

where S is a positive infinitesimal. For the clean surface, the eigenvalues Ej form a continuum 
of constant density of states po = 0.045 states/(eV site), using an effective mass 0.38 times 
the electron mass m e |0|. For a hard wall elliptic corral, the Ej are discrete. For a finite 
confinement potential, the eigenstates inside the corral become resonances and the spectrum 
is continuous again. However, from results borrowed mainly from nuclear physics, we know 
that G° inside the corral can be expressed as a discrete sum of complex poles |^|. Thus, the 
form of Gc for the hard wall corral is retained, but 5 becomes finite. Here we retain the values 
of Ej and <fj(r) of the hard wall corral calculated before pj up to an energy of 1 eV above 
e F - This approximation is justified by comparison with the Green function which results for 
a scattering potential consisting of delta functions regularly spaced at the boundary of the 
ellipse [B]. For all continuum states around a resonance, the wave function inside the corral 
is essentially the same, differing only outside the corral. In adition, with an effective mass 
0.378 m e , the hard wall results for Ej around e F reproduce the experimental results. For 
simplicity, we take 5 independent of j. 




Figure 1: Change of the spectral density of the mixed state / for (a) clean surface, (b-d) 
elliptical corral with impurity at the left focus, and the indicated values of <5 and q. The 
tip is at the impurity position except for the dotted lines which correspond to the empty 
focus in the corral 



The impurity Green function is given by: 

G- d 1 {uj)=Lu-E d ~V 2 G < l{R l ,R l ,Lj)-Y 1 {Lo), (5) 

where in our approach, E d is a renormalized d level, and £(o>) is the self energy calculated 
in second order in U M. Because of the limitations of the perturbative approach, we have 
chosen a moderate value of U = 1 eV, and E d — tF (near the symmetric case). We have 
taken V = 0.64 eV, to lead to a half width at half maximum near Tk = 4 meV for the 
spectral density pd(w) of the impurity on a clean surface Q. The resulting impurity resonant 
level width is A ~ 58 meV. For this value of U/A, the position of the peak is not adequately 
given by perturbation theory. Then we used independent values of E d for the impurity on 
the clean surface or inside the corral, to control this position. The other parameters of the 
model are the same in both situations. The width 8 is the only free parameter to control the 
line shape for the impurity inside the corral. 

Results In Fig. 1 (a), we show Apf, proportional to AdI/dV, at the impurity position for 
the clean surface and several values of q. The result for q = is symmetrical and proportional 
to pd{w). The experimental AdI/dV is larger for positive sample bias, corresponding to 
positive energy u). This asymmetry can be achieved increasing q. However, the experimental 
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line shape shows a relative maximum also at negative u), in contrast to our result. This 
feature seems difficult to explain. Except for this fact, the observed line shape is qualitatively 
reproduced for q between 0.03 and 0.04. 

In Fig. I (b), (c) and (d) we show Apf for an elliptical corral with eccentricity e = 1/2, 
and semimajor axis a = 71.3 A §, with an impurity placed at the left focus, for several 
values of 8 and q, and for both foci. If 8 < Tk, the Kondo resonance at €f is absent. Instead, 
both pd and p/ present two narrow peaks at both sides of the Fermi energy. In the absence 
of the impurity, pf has a peak at which corresponds to the state 42 {ep = £42). As 
a consequence, the difference Apf present two narrow peaks above zero and a depression 
with negative Apf (see Fig. 1 (b)). For 8 3> Tk (Fig. 1 (c) and (d)), a Kondo resonance 
quite similar to that of the clean surface is formed. However, due to the asymmetry of the 
amplitude of the wave functions of the corral above and below €f at the impurity site, the line 
shape is asymmetric for q = 0, being higher at negative energies. This asymmetry is opposite 
to the experimentally observed in the case of the clean surface. Thus, while a positive value 
of q renders the line shape asymmetric for a clean surface, it restores a rather symmetric line 
shape in the case of the corral, as experimentally observed. 

As 8 increases, the magnitude of the depression observed at the empty focus decreases. 
This is due to an increasing negative interference of the corral state 42, which is even under 
interchange of both foci, with other odd states, whose weight at cf increase with 8 M. For 
8 ~ 40 meV, the main aspects of the observed space and voltage dependence of AdI /dV are 
reproduced. 



Conclusion An Anderson model in which the magnetic impurity is hybridized with surface 
states, which are extended in the case of a clean surface or resonances in presence of a quantum 
corral, is able to reproduce the qualitative features of the observed change in the voltage 
dependence of the differential conductance AdI /dV (in particular its asymmetry) in both 
cases. This requires a small direct tunneling between the impurity and probing tip, when 
the tip is above the impurity. The rather symmetric line shape observed above the impurity 
at the focus in the quantum corral is the result of the compensation of two effects which 
alone would produce an asymmetry in opposite directions: the direct tunneling between the 
impurity and tip, and the particular electronic structure of the resonances inside the corral, 
above and below the Fermi energy. Work to study the line shape in other physical situations, 
for example for a mirage out of focus [^, ^, |l(J is in progress. 
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